Electron-acoustic solitary waves in the presence of a suprathermal electron component 
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The nonlinear dynamics of electron-acoustic localized structures in a collisionless and unmagne- 
tized plasma consisting of "cool" inertial electrons, "hot" electrons having a kappa distribution, and 
stationary ions is studied. The inertialess hot electron distribution thus has a long-tailed suprather- 
mal (non-Maxwellian) form. A dispersion relation is derived for linear electron-acoustic waves. They 
show a strong dependence of the charge screening mechanism on excess suprathermality (through 
k) . A nonlinear pseudopotential technique is employed to investigate the occurrence of stationary- 
profile solitary waves, focusing on how their characteristics depend on the spectral index k, and 
the hot-to-cool electron temperature and density ratios. Only negative polarity solitary waves are 
found to exist, in a parameter region which becomes narrower as deviation from the Maxwellian 
(suprathermality) increases, while the soliton amplitude at fixed soliton speed increases. However, 
for a constant value of the true Mach number, the amplitude decreases for decreasing k. 

PACS numbers: 52.35.Sb, 52.35.Mw, 52.35.Fp, 72.30.+q 



I. INTRODUCTION 

Electron-acoustic waves may occur in plasmas charac- 
terized by a co-existence of two distinct electron popu- 
lations, here referred to as "cool" and "hot" electrons. 
These are electrostatic waves of high frequency (in com- 
parison with the ion plasma frequency), propagating at 
a phase speed which lies between the hot and cool elec- 
tron thermal velocities. On such a fast (high frequency) 
dynamical scale, the positive ions may safely be assumed 
to form a uniform stationary charge background simply 
providing charge neutrality, yet playing no essential role 
in the dynamics. The cool electrons provide the inertia 
necessary to maintain the electrostatic oscillations, while 
the restoring force comes from the hot electron pressure. 

A matter of importance in electrostatic wave propaga- 
tion (although inevitably overlooked in fluid plasma mod- 
els) is Landau damping, which becomes stronger when 
the phase velocity approaches the thermal velocity of ei- 
ther electron component, thus the wave can propagate 
in the plasma only within a restricted range of param- 
eter values. It turns out that electron-acoustic waves 
are weakly damped for a temperature ratio Tc/Th < 0.1 
and provided that the cool electrons represent an in- 
termediate fraction of the total electron density: 0.2 < 
"n-c/inc + nfi) ^ 0.8. [iHJI The wavenumber k to minimize 
damping lies roughly between 0.2A^^ and 0.6A^J. (where 
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Xdc is the cool electron Debye length). These results 
on bi-Maxwellian plasmas were later extended to include 
the effect of the excess suprathermality of the hot elec- 
trons p\ (a physical feature to be discussed below). It 
was found that excess suprathermal electrons do cause a 
modification of the damping curves, but the overall qual- 
itative conclusion remains unchanged: electron-acoustic 
waves survive Landau damping over a wide range of pa- 
rameter values. (5| However, care must be taken in the 
choice of plasma configuration when studying nonlinear 
electron-acoustic structures, so as to ensure that one is 
considering a region of parameter space in which Landau 
damping is minimized. 

Electron-acoustic waves occur in laboratory experi- 
ments M H 8'iid space plasmas, e.g., in the Earth's bow 
shock 0-[ip| and in the auroral magnetosphere [l|, lll| . 
They are associated with Broadband Electrostatic Noise 
(BEN), a common high-frequency background activity, 
regularly observed by satellite missions in the plasma 
sheet boundary layer (PSBL) |12h14| . BEN emission in- 
cludes a series of isolated bipolar pulses, within a fre- 
quency range from ~ 10 Hz up to the local electron 
plasma frequency (^ 10 kHz) [l^- This clearly suggests 
that BEN is related to electron dynamics rather than to 
the ions 12, Jl4] . 

In the standard bi-Maxwellian picture, the two electron 
species would each be assumed to be in a (different) ther- 
mal Maxwellian distribution, parameterized via two dis- 
tinct temperature values, Tc and T^, respectively |15l - [l7| . 
Contrary to this picture, space and laboratory plasmas 
often possess an excess population of suprathermal elec- 
trons, a fact which is reflected in a power law distribution 
at high velocity (above the electron thermal speed). This 
excess suprathermality phenomenon is well modeled by 
a generalized Lorentzian or K-distribution |18l - [23 |. The 
common form of the isotropic (three-dimensional) gener- 



alized Lorentzian or K-distribution function is given by 
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where no is the equihbrium number density of the elec- 
trons, V the velocity variable, and 9 the most proba- 
ble speed, which acts as a characteristic "modified ther- 
mal speed", and is related to the usual thermal speed 

vth,e = {2kBTe/mef''^ by 6* = vth,e [[k - 1)1 n] ^'^. Here 
Ub is the Boltzmann constant, rUe the electron mass and 
Te the temperature of an equivalent Maxwellian having 
the same energy content. J22| The term involving the 
Gamma function (F) arises from the normalization of 
/k(w), viz., J ff^{v)(fiv = riQ. Here, suprathermality is 
denoted by the spectral index k, with k > |, for real- 
ity of the most probable speed, 9. [23| Low values of n 
are associated with a significant number of suprathermal 
particles; on the other hand, for k — > oo a Maxwellian 
distribution is recovered. 

The K-distribution was first applied to model veloc- 
ity distributions observed in space plasmas that were 
Maxwellian-like at lower velocities, but had a power-law 
form at higher speeds [23| , and was later applied in a va- 
riety of studies, successfully fitting many real space ob- 
servations, e.g., 0, [21I, [2J, [23 ■ Typical k values usually 
lie in the range 2 < k < 6. For example, observations 
in the earth's foreshock satisfy 3 < Kg < 6, [0] mea- 
surements of plasma sheet electron and ion distributions 
yield k^ = 4.7 and Kg = 5.5 (here, e denotes electrons 
and i ions), [2J| and coronal electrons in the solar wind 
are modeled with 2 < Kg < 6 [2^. Recent observations 
of the radial distribution of the electron population in 
Saturn's magnetosphere also point towards a kappa dis- 
tribution (we — 2.9 - 4.2) [26|. Therefore, we focus our 
interest in the following on the range 2 < k < 6; in fact, 
the Maxwellian limit is already practically attained for 
values above k ~ 10. 

A linear analysis of electron- acoustic waves was first 
carried out by assuming an unmagnetized Maxwellian 
homogeneous plasma, which exhibited a heavily damped 
acoustic-like solution in addition to Langmuir waves and 
ion-acoustic waves. [23] Those early results were later ex- 
tended to take into account the effect of excess suprather- 
mal particles [a, [28| , whose presence in fact results in an 
increase in the Landau damping at small wavenumbers, 
in particular when the hot electron component is dom- 
inant [5|, |29|. Studies of linear and nonlinear electron- 
acoustic waves in plasmas with nonthermal electrons have 
received a great deal of interest in recent years [30l - t34| . 
Negative potential solitary structures were shown to ex- 
ist in a two-electron plasma, either for Maxwellian [30| 
or for nonthermal [3l|,[3a] hot electrons. Interestingly, ei- 
ther incorporation of finite inertia |32l . |33| or the addition 
of a beam component [3a [33| may lead to the existence 
of positive and negative potential solitons. A recent in- 
vestigation has established the properties of modulated 



electron-acoustic wavepackets in kappa-distributed plas- 
mas, and has studied the effect of su pra thermality on the 
amplitude (modulational) stability, [sl] 

In this paper, we study the linear and nonlinear dy- 
namics of electron-acoustic waves in a plasma consisting 
of cool adiabatic electron and hot ^-distributed electrons, 
in addition to stationary ions. The paper is organised as 
follows. An electron-plasma-fluid model is presented in 
Section ini In Section Hill a linear dispersion relation is 
derived and discussed. In Section |IV1 the Sagdeev pseu- 
dopotential method is employed to investigate the occur- 
rence of stationary profile electrostatic solitary waves. In 
Section|Vl we depict the existence domain of the electron- 
acoustic solitary waves. Section |VT] is devoted to a para- 
metric investigation of the form of the Sagdeev pseudopo- 
tential and of the characteristics of electron-acoustic soli- 
tary waves. Our results are summarized in Section [VIII 



II. MODEL EQUATIONS 

We consider a plasma consisting of three components, 
namely a cool electron-fluid (at temperature Tc ^ 0), an 
inertialess hot electron component with a nonthermal (k) 
velocity distribution, and uniformly distributed station- 
ary ions. 

The cool electron behavior is governed by the continu- 
ity equation. 
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The pressure of the cool electrons is governed by 
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Here ric, u^ and pc are the number density, the velocity 
and the pressure of the cool electron fluid, is the elec- 
trostatic wave potential, e the elementary charge, and 
J — if + 2)// denotes the specific heat ratio (for / de- 
grees of freedom). We shall assume 7 = 8 (viz., / = 1 in 
ID) for the adiabatic cool electrons. 

We assume the ions to be stationary (immobile), i.e., 
in a uniform state rii = rii.o — const, (where n^.o is the 
undisturbed ion density) at all times. In order to obtain 
an expression for the number density of the hot electrons, 
rih , based on the n distribution ([T]) , one may integrate Eq. 
(HI) over the velocity space, to obtain 20] 



nh{(t)) = Uhfi 1 - 



eip 



ksThin- 2) 



-1/2 



(5) 



where Uh.o and Th are the equilibrium number density 
and "temperature" of the hot electrons, respectively, and 



K is the spectral index measuring the deviation from ther- 
mal equilibrimxi. 

The densities of the (K-distributed) hot electrons, the 
adiabatic cool electrons, and the stationary ions are cou- 
pled via Poisson's equation: 
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where Eq is the permittivity constant, Uh and rii are the 
number density of hot electrons and ions, respectively. 
At equilibrium, the plasma is quasi-neutral, so that 



rico + ntfi = Zni,Q , 



(7) 



implying Znifi/ricfi = 1-1-/3, where we have defined the 
hot-to-cool electron density ratio 
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According to Ref. X Landau damping is minimized 
in the range 0.2 < rico/incfi + nhfi) < 0.8, implying 
0.25 ^ (3 ^ 4. This is our region of interest in what 
follows, as nonlinear structures will not be sustainable 
for plasma configurations for which the linear waves are 
strongly damped. 

Scaling by appropriate quantities, we obtain the nor- 
malized set of equations 
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Here, n, u and p denote the cool electron fluid den- 
sity, velocity and pressure variables normalized with re- 

spect to nc,o, Cth = [fesTl/me] and nc^oksTc, respec- 
tively. Time and space were scaled by the plasma period 
'^pc = ("-c,oe^/£o"T-e)~^^^ and the characteristic length 
Ao = {eokBTh/ncfie^y^^ , respectively. Finally, is the 
wave potential scaled by ksTh/e. We have defined the 
temperature ratio of the cool to the hot electrons as 



a = Tc/Th. 



III. LINEAR WAVES 



(13) 



As a first step, we linearize Eqs. (^t- P^ . to study 
small-amplitude harmonic waves of frequency uj and 
wavenumber k. The linear dispersion relation for 
electron-acoustic waves then reads: 
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where a/So- is essentially the (normalized) cool electron 
thermal velocity. After taking account of differences in 
normalization, this agrees with the form found in [39|. 

We note the appearance of a normalized K-dependent 
screening factor (scaled Debye wavenumber) fc^i^K in the 
denominator, defined by 
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Since this is the inverse of the (Debye) screening length, 
we notice that the latter in fact decreases due to an ex- 
cess in suprathermal electrons (i.e., Xd,k < Xd,oo for any 
finite value of k). This is in agreement with Refs. |40| - |43 : 
note also the discussion in Ref. ISSl 

From Eq. ([T4)) . we see that the frequency a;(/c), and 
hence also the phase speed, increases with higher tem- 
perature ratio a ~ Tc/Th- However, this is usually a 
small correction to the dominant first term on the right- 
hand side of (|15l) . For large wavelength values (small 
k ^ kD,K), the phase speed is given by 
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while on the other hand, the thermal contribution is dom- 
inant for high wavenumber k ^ fcc.Kj i-e. 



Vph ~ (3cr) 



1/2 



(17) 



However, we should recall from kinetic theory [iHS Q 
that both for very long wavelengths {kXoc '^ 0.2) and 
very short wavelengths {kXoc > 0.6), the wave is strongly 
damped, and thus these limits may be of academic inter- 
est only. The mode is weakly damped only for interme- 
diate wavelength values, where its acoustic nature is not 
manifest. [IHl H, Hi Here, Xdc = {eokBTc/ucfie^y/^ is 
the cool electron Debye length. 

Restoring dimensions for a moment, the dispersion re- 
lation becomes 
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where ctc = {kBTc/nieY^'^ is the cool electron thermal 
speed and Xoh is the hot electron Debye length defined 
by 
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It appears appropriate to compare the above results 
with earlier results, in the linear regime. First of all, 
we note that Ref. y has adopted a kinetic description 
of electron- acoustic waves in suprathermal plasmas. For 
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IV. NONLINEAR ANALYSIS FOR LARGE 
AMPLITUDE SOLITARY WAVES 

Anticipating constant profile solutions, we shall con- 
sider Eqs. ([^- ([T^ in a stationary frame traveling at a 
constant normalized velocity M (to be referred to as the 
Mach number), implying the transformation ^ = x — Mt. 
The space and time derivatives are thus replaced by 
djdx = d/dS, and d/dt — ~Md/d£^, respectively, so Eqs. 
©-([Ill) take the form: 
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FIG. 1. (Color online) Dispersion curve for harmonic (linear) 
electron-acoustic waves. Upper panel: The variation of the 
dispersion curve for different values of k is depicted. Curves 
from top to bottom: k = 10 (solid), 4 (dashed), 3 (dot- 
dashed), and 2 (dotted curve). Here, a — 0.01 and P — 1. 
Bottom panel: Variation of the dispersion curve for different 
values of fi. Curves from top to bottom: (5 = 0.5 (solid), 
1 (dashed), and 3 (dot-dashed curve). Here, (j — 0.01 and 
K = 3. 



this purpose, Eq. (|18p may be cast in the form. 
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Here Cl. 






("c.o/«m)[('«-3/2)/(k-1/2)]c; 



is the electron acoustic speed in a suprathermal plasma, 
while the effective shielding length is given by A^^ = 
[€okBTh/nh,oe'^][{K.-i/2)/{n-l/2)]. Our Eq. ([101 above 
agrees with relation (3) in Ref. i5! (upon considering the 
limit \dc — >■ therein). In the limit k — >■ oo (Maxwellian 
distribution), Eq. (|20l) recovers precisely Eq. (1) in Ref. 
l36l . upon considering the limit Adc — > 0. Furthermore, 
one recovers exactly Eq. (5) in Ref. \\A for Maxwellian 
plasma in the cold-electron limit {Tc= 0). 

In Figure [1] we depict the dispersion curve of the 
electron-acoustic mode, showing the effect of varying the 
values of the spectral index k and the density ratio /3. It 
is confirmed numerically that the phase speed (uj/k) in- 
creases weakly with a reduction in suprathermal particle 
excess, as the Maxwellian is approached, and that there 
is a significant reduction in phase speed as the plasma 
model changes from one in which the cool electrons dom- 
inate, to one which is dominated by the hot electron den- 
sity. 
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We assume that the equilibrium state is reached at both 
infinities (^ — >■ ±oo). Accordingly, we integrate and apply 
the boundary conditions n = l,p=l,u = and (/> = 
at ±oo. One thus obtains 
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u = M - (mV2(/) - 3nV + 3cr) ' , (26) 



and 
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Combining Eqs. (|25p -(|27 |) . we obtain the following 
biquadratic equation for the cool electron density. 



The solution of Eq. (|28l) may be written as 
1 
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From the boundary conditions, n = 1 at (j> ~ 0, it follows 
that the negative sign must be taken in Eq. (P^)) . Fur- 
thermore, we shall assume that M > a/3ct, i.e., that the 
cool electrons are supersonic, while the hot electrons are 
subsonic, thus we require that M < 1. 

Reality of the density variable imposes the requirement 
2(f> + (M— \/3ct) > 0, which implies a limit on the elec- 
trostatic potential value |(/)max| — | (M— -s/Sct) associ- 
ated with negative solitary structures (positive electric 
potentials, should they exist, satisfy the latter condition 
automatically, and are thus not limited). 

Substituting the density expression (|29p - (PT|) into 
Poisson's equation ([M]) and integrating, yields the 
pseudo-energy balance equation for a unit mass in a con- 
servative force field, if one defines ^ as "time" and </> as 
"position" variable: 
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where the Sagdeev pseudopotential ^((/)) is given by 
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V. SOLITON EXISTENCE DOMAIN 

We next investigate the conditions for existence of soli- 
tons. First, we need to ensure that the origin aX (j) — 
is a root and a local maximum of ^ in Eq. (I33p. i.e., 
*(0) = 0, *'(0) = and *"(0) < at </) = [S-lil, 
where primes denote derivatives with respect to (f>- It is 
easily seen that the first two constraints are satisfied. We 
thus impose the condition 
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Eq. (j34p provides the minimum value for the Mach num- 
ber, Ml, i.e. : 
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Clearly, Mi is the (normalized) electron-acoustic phase 
speed - cf. Eq. (|16l) . It is thus also related to Debye 
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FIG. 2. (Color online) Variation of the lower limit Mi (lower 
curves) and the upper limit M2 (upper curves) with the hot- 
to-cold electron density ratio /3 for different values of the tem- 
perature ratio a. SoUtons may exist for values of the Mach 
number M in the region between the lower and the upper 
curve(s) of the same style/color. Curves: (a-b) a = 0.01 
(solid), 0.02 (dashed), and 0.04 (dot-dashed), and (c) a = 0.01 
(solid), 0.1 (dashed), and 0.2 (dot-dashed) Here, we have 
taken: (a) k — 2, (b) k — 100 (quasi-Maxwellian), and (c) 
K^ 1.65. 



screening via the screening parameter A^i.k in (|15p . asso- 
ciated with the hot K-distributed electrons. We deduce 
that soliton solutions are super-acoustic. For Maxwellian 
hot electrons {k — >■ 00) and cold "cool" electrons {a — 0), 
we obtain Mi = l//?^/^, thus recovering the normalized 
phase speed for electron-acoustic waves in a Maxwellian 
plasma. [5| The lower Mach number limit. Mi, increases 
with Tc (via a), and decreases for lower values of k (large 
excess of suprathermal electrons), and hence the sound 
speed in suprathermal plasmas is reduced, in comparison 
with Maxwellian plasmas (k — > 00). 

An upper limit for M is found through the fact that 
the cool electron density becomes complex aX (j) = <t>max, 
and hence the largest soliton amplitude satisfies F2 (Af ) = 
^(0)|</,=0„,.,^ > 0. This yields the following equation for 




FIG. 3. (Color online) Variation of the lower limit Mi (lower 
curves) and the upper limit M2 (upper curves) with the 
suprathermality parameter k for different values of the tem- 
perature ratio a (upper panel), and density ratio j3 (bottom 
panel). Solitons may exist for values of the Mach number 
M in the region between the lower and upper curves of the 
same style/color. Upper panel: a = 0.01 (solid curve), 0.02 
(dashed), and 0.04 (dot-dashed). Here, we have taken P = 2. 
Lower panel: /3 = 1.8 (solid), 2 (dashed), 2.2 (dot-dashed), 
and 5 (solid circles). Here, a = 0.02. 
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FIG. 4. (Color online) The pseudopotential '!'((?!') (upper 
panel) and the associated solution (electric potential pulse) 
(\> (lower panel) are depicted versus position ^, for differ- 
ent values of the temperature ratio a . We have taken: 
a = 0.01 (solid curve), 0.02 (dashed curve), and 0.04 (dot- 
dashed curve). The other parameter values are /3 = 1.3, 
K = 2.5 and M = 0.75. 



the upper limit in M: 
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Solving Eq. (jJH) provides the upper limit M2(k, /3, cr) 
for acceptable values of the Mach number for solitons to 
exist. 

For comparison, for a Maxwellian distribution (here 
recovered as k — >■ cx)), the constraints reduce to 
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The latter equation provides the upper limit Af2, while 
the lower limit becomes A/i = (1//3 -t- SaY^^. 



In the opposite limit of ultrastrong suprathermality, 
i.e., K — >■ 3/2, the Mach number threshold approaches a 
non-zero limit Afi = \/3ct, which is essentially the ther- 
mal speed, as noted above (recall that M > ^/Sa by 
assumption). The upper limit M2 is then given by 



F2(Af) = -^(1 + (3) (a/Z-n/S^) + M^+a 



-A/3/2 (3(7)^/" = 0. 
3 ^ ' 



(39) 



Interestingly, the two limits A/i and Af2 both tend to 
the same limit as k — > 3/2, namely, v3o^, where the soli- 
ton existence region vanishes, as the kappa distribution 
breaks down. 

We have studied the existence domain of electron- 
acoustic solitary waves for different values of the param- 
eters. The results are depicted in Figs. [2H3l Solitary 
structures of the electrostatic potential may occur in the 
range A/i < M < A/2, which depends on the parameters 
/?, K, and a. We recall that we have also assumed that 
cool electrons are supersonic (in the sense M > \/3ct) 
|43l - l45{ , and the hot electrons subsonic {M < 1), and 
care must be taken not to go beyond the limits of the 
plasma model. 

The interval [A/i,Af2] where solitons may exist is de- 
picted in Fig. [21 in two opposite cases: in (a) and (c) two 



very low, and in (b) one very high value of k. We thus 
see that for both a quasi-Maxwellian distribution and one 
with a large excess suprathermal component of hot elec- 
trons, both Ml and M2 decrease with an increase in the 
relative density parameter /3 for fixed k and soliton speed 
M . Further, the upper limit falls off more rapidly, and 
thus the existence domain in Mach number becomes nar- 
rower for higher values of the hot-to-cool electron density 
ratio. Comparing the two frames (a) and (b) in Fig. [2l we 
immediately notice that suprathermality (low k) results 
in solitons propagating at lower Mach number values, a 
trend which is also seen in Fig. [5};. Another trend that 
is visible in Figs. [2H3^ is that increased thermal pressure 
effects of the cool electrons, manifested through increas- 
ing a, also lead to a narrowing of the Mach number range 
that can support solitons. Finally, we note that for /3 ~ 1, 
the upper limit found from Eq. (j36l) rises above the limit 
M = 1 required by the assumptions of the model, and 
the latter then forms the upper limit. 

Interestingly, in Figs. EHSlthe existence region appears 
to shrink down to nil, as the curves approach each other 
for high /3 values. This is particularly visible in Fig. [2t, 
for a very low value oi k {k — 1.65). This is not an un- 
expected result, as high values of /? are equivalent to a 
reduction in cool electron relative density, which leads to 
our model breaking down if the inertial electrons vanish. 
We recall that a value /? > 4 is a rather abstract case, 
as it corresponds to a forbidden regime, since Landau 
damping will prevent electron-acoustic oscillations from 
propagating. Similarly, a high value of the temperature 
ratio, such as ct = 0.2, takes us outside the physically 
reasonable domain. Nevertheless, as it appears that the 
lower and upper limits in M approach each other asymp- 
totically for high values of /3, we have carried out calcu- 
lations for increasing /3, up to /? = 100 for cr = 0.2 as an 
academic exercise, and can confirm that the two limits 
do not actually intersect. 

Figure |3] shows the range of allowed Mach numbers as a 
function of k, for various values of the temperature ratio 
a. As discussed above, increasing n towards a Maxwellian 
distribution (k — >■ c») broadens the Mach number range 
and yields higher values of Mach number. On the other 
hand, both upper and lower limits decrease as the limit- 
ing value K — >■ 3/2 is approached. The qualitative conclu- 
sion is analogous to the trend in Fig. [2] stronger excess 
suprathermality leads to solitons occurring in narrower 
ranges of M. Furthermore, as illustrated in Figs. [5] and 
[5^, the Mach number threshold Mi approaches the upper 
limit M2 for high values of a and (3: both increased hot- 
electron density and cool-electron thermal effects shrink 
the permitted soliton existence region. 

Figure [3}d depicts the range of allowed Mach numbers 
as a function of n for various values of the density pa- 
rameter /3 (for a fixed indicative a value). We note that 
both curves decrease with an increase in /3. Although it 
lies in the damped region, we have also depicted a high 
/3 regime for comparison (solid-crosses curve). 

We conclude this section with a brief comparison of our 
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FIG. 5. (Color online) (a) The pseudopotential *I'('^) and the 
associated solutions: (b) electric potential pulse </>, (c) density 
n, and (d) velocity u are depicted versus position S,, for differ- 
ent K. We have taken: k — 2.5 (solid curve), 3 (dashed curve), 
and 3.5 (dot-dashed curve). The other parameter values are: 
a = 0.02, /3 = 1.6, and M = 0.75. 



work with that of Ref. [39|. The latter did not consider 
existence domains at all, let alone their dependence on 
plasma parameters, but merely plotted some Sagdeev po- 
tentials and associated soliton potential profiles for cho- 
sen values of some of the parameters, so as to extract 
some trends. En passant, there is indirect mention of an 
upper limit in M, in that it is commented that as increas- 
ing values of M are considered, at some stage solitary 
waves cease to exist. [33 
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FIG. 6. (Color online) The dependence of the pulse ampli- 
tude l^ml on the Mach number-to-sound-speed ratio M/M\ 
is depicted, for different values of k. From top to bottom: 
K. — 100 (solid curve); 10 (dashed curve); 4 (dot-dashed 
curve); 3 (crosses); 2.5 (solid circles); 2 (solid squares). Here, 
a = 0.01 and /3 = 1.3. 
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FIG. 7. (Color online) (Upper panel) The pseudopotential 
'3/(0) vs. (\> and (lower panel) the associated electric potential 
pulse (^ vs. 5 are depicted, for different values of the hot-to- 
cold electron density ratio fi. From top to bottom: /? = 1.3 
(solid curve) ; 1.5 (dashed curve) ; 1.7 (dot-dashed curve) . Here 
a = 0.01, K = 2.5 and M = 0.75. 



VI. SOLITON CHARACTERISTICS 

Having explored the existence domains of electron- 
acoustic solitons, subject to the constraints of the plasma 
model, we now turn to consider aspects of the soliton 
characteristics. We have numerically solved Eq. (15^ for 



different representative parameter values, in order to in- 
vestigate their effects on the soliton characteristics. We 
point out that, varying different parameters, we have 
found only negative potential solitons, regardless of the 
value of K considered. This is not altogether unexpected, 
as it has been found in a number of examples that, in 
contrast to the Cairns model, [46] , the kappa distribution 
does not lead to reverse polarity acoustic solitons. i2Qi,i47| 

Figure |4] shows the variation of the Sagdeev pseudopo- 
tential ^(0) with the normalized potential 0, along with 
the associated pulse solutions (soliton profiles) , for differ- 
ent values of the temperature ratio, a = T^/Th (keeping 
/3 = 1.3, K, = 2.5 and Mach number M = 0.75, all fixed). 
The Sagdeev potential well becomes deeper and wider as 
a is increased. We thus find associated increases in the 
soliton amplitude and in profile steepness (see Fig. |3b)- 
Thermal effects therefore are seen to amplify significantly 
the electric potential disturbance at fixed M . 

Figure [S^ shows the Sagdeev pseudopotential *((/i) for 
different values of n. The electrostatic pulse (soliton) 
solution depicted in Fig. [5}d is obtained via numerical in- 
tegration. The pulse amplitude \4>m\ increases for lower 
K, implying an amplification of the electric potential dis- 
turbance as one departs from the Maxwellian. Once the 
electric potential has been obtained numerically, the cool- 
electron fluid density (Fig. [St) and velocity disturbances 
(Fig. [SJi) are determined algebraically. Both these dis- 
turbances are positive in this case, and again, for lower 
K values, the profiles reflecting the compression and the 
increase in velocity are steeper but narrower. 

We recall that as various parameter values are varied, 
the true acoustic speed in the plasma configuration. Mi , 
also varies. As solitons are inherently super-acoustic, it 
is clear that the effect of a changing true acoustic speed 
could mask other dependences. Hence it is also desir- 
able to explore soliton characteristics as a function of 
the propagation speed M, measured relative to the true 
acoustic speed. Mi. This ratio, M/Mi, thus represents 
the "true" Mach number. It has been shown that for any 
plasma made up of barotropic fluids, arbitrary amplitude 
solitons satisfy d'^/dM < 0, [isl-ls^ from which it follows 
that d(j)rn/dM > 0, where (j)m is the soliton amplitude. 
Thus one expects that the soliton amplitude is an in- 
creasing function of M/Mi. This is true for both KdV 
solitons (small amplitudes, propagating near the sound 
speed) - "taller is faster"- and, in principle, also for fully 
nonlinear (Sagdeev) pulses (where the soliton character- 
istics can only be found numerically. [43, |51|) In Fig. [51 
we have plotted the soliton amplitude \4>ui\ as a function 
of the ratio M/Mi , for a range of values of the parameter 
K. Clearly, the amplitude increases linearly with M /Mi 
for all values of k. The two plots for k < 2.5 both cover 
the full range up to M — M2- However, although we 
deduce from earlier figures that M2 increases with k, we 
see that the endpoints of the plots for k > 3 occur at 
decreasing values of MjM\ , and indeed decreasing max- 
imum amplitudes (\)m- That occurs as, for the chosen 
values of /? and cr, M2 exceeds unity for k > 3, and we 



have truncated the curves at the point where M = 1, to 
remain within the range defined by the plasma modeL 

The effect of the hot-to-cool electron density ratio, /3 
on the soliton characteristics is shown in Fig. [71 We 
see that the soliton excitations are amplified and profiles 
steepened (the Sagdeev potential well becomes wider and 
deeper) , as the density of the hot (nonthermal) electrons 
is increased (i.e., for higher /3), viz., keeping k, a and M 
fixed. Furthermore, an increase in the number density 
of the hot electrons also leads to an increase in the per- 
turbation of both density n, and velocity u of the cool 
electrons (figure omitted). 



VII. CONCLUSION 

In this article, we have performed a thorough linear 
and nonlinear analysis, from first principles, of electron 
acoustic excitations occurring in a nonthermal plasma 
consisting of hot K-distributed electrons, adiabatic cool 
electrons, and immobile ions. 

First, we have derived a linear dispersion relation, and 
investigated the dependence of the dispersion character- 
istics on the plasma environment (degree of 'suprather- 
mality' through the parameter k, plasma composition, 
and thermal effects). 

Then, we have employed the Sagdeev pseudopotential 
method to investigate large amplitude localized nonlin- 
ear electrostatic structures (solitary waves), and to deter- 
mine the region in parameter space where stationary pro- 
file solutions may exist. Only negative potential solitons 
were found. The existence domain for solitons was shown 
to become narrower in the range of solitary wave speed, 
with an increase in the excess of suprathermal electrons in 
the hot electron distribution (stronger 'suprathermality', 
lower K value) . The dependence of the soliton character- 
istics on the hot electron number density (through the 
parameter /3) and on the hot-to-cool electron tempera- 
ture ratio a, were also studied. A series of appropriate 
examples of pseudopotential curves and soliton profiles 
were computed numerically, in order to confirm the pre- 
dictions arising from the study of existence domains. 

It may be added that ionic motion/inertia, here ne- 
glected, may also be included for a more accurate de- 
scription, but is likely to have only minor quantitative 
effects. 

We note, for completeness, that very recently two re- 
lated papers have appeared with a scope apparently sim- 
ilar to that of the present article, viz., Refs. [3^, [S^]- A 
word of comparison may therefore be appropriate here, 
for clarity. The latter authors 52] indicate that they 
are using a form of kappa distribution from one of the 
pioneering papers in the field. [53] Unfortunately, their 
expression for the characteristic speed 9 does not agree 
with the standard expression [Sj], and thus the hot elec- 
tron density does not take the usual form, Eq. (5). |20| 
Further, they use values of k > 0.6, i.e., well below the 
standard K-distribution cut-off of 3/2. Hence their results 



do not apply to the standard form of k distribution. [22] 

As regards the other paper, [33] it did not consider 
existence domains at all (apart from an indirect mention 
of an upper limit in M, as commented on in Section [Vj 
above). Further, no account is taken in the paper of 
the possible effects of Landau damping on sustainable 
nonlinear structures, and a number of the figures relate 
to values of /3 (called a in the paper) which lie in the 
unphysical, damped range. 

We note that Ref. [33 has also carried out a "small 
amplitude" calculation yielding double layers. However, 
when evaluated numerically, these turn out to be well be- 
yond the range of small amplitude, and that raises some 
doubts about the validity of the results (their Figures 8- 
10). Further, let us take together their Figures 2 and 8, 
and consider the case of k = 3 and a — 0.2 (i.e., our 
/3). The latter is, of course, a value for which the linear 
wave is likely to be strongly Landau damped. It ap- 
pears from the figures that a soliton occurs at M = 1.1 
with amplitude ~ 0.9, while a double layer occurs for 
M — 2.0 with amplitude ~ 0.7. This combination of data 
does not satisfy the analytically-proven requirement that 
d'^/dM < 0, [isl - lsoj . There is no obvious reason why 
that should be the case, and there thus seems to be an 
error in at least one of these two figures. 

Finally, we point out that we have not sought double 
layers in our calculations. However, we would be sur- 
prised if they did occur, as they are usually found as the 
upper limit to a sequence of solitons for a polarity for 
which there is no other limit. In this case there clearly is 
an upper limit for negative solitary waves, arising from 
the constraint F2{M) = 0. It is in principle possible for 
a double layer to occur at a lower value of Af , and be 
followed by larger amplitude solitons at higher A/, until 
the upper fiuid cutoff such as a sonic point or an in- 
finite compression cutoff is reached. |54J ] However, such 
behaviour depends on the the Sagdeev potential having a 
fairly complicated shape, with subsidiary local maxima, 
and we have not observed these for this model. 

We are not aware of any experimental studies with 
which these theoretical results may be directly compared. 
However, it has previously been shown that wave data 
may be used to obtain an estimate for k, thus acting as 
a diagnostic for the distribution function [13, [S^] . 

Similarly, in this case there are a number of indica- 
tors amongst our results which experimenters may wish 
to consider when interpreting observations. Thus, for 
instance, a lower normalized phase velocity of the lin- 
ear electron-acoustic wave than would be predicted by a 
Maxwellian model (see Fig. [T]) could be used to evaluate 

K. 

Secondly, from Fig. [5] one sees that in 1ow-k plas- 
mas the range of normalized soliton speeds is both nar- 
rower and of larger value than one would expect for a 
Maxwellian. Thus, if solitons are found with normalized 
speeds around M ~ 0.4, these can be understood only by 
allowing for additional suprathermal electrons (lower n) . 
Further, from Fig. [5] it follows that Maxwellian electrons 
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give rise to a cutoff in the density ratio /3 — 1 , and hence 
sohtons observed in such plasmas can only be explained 
in terms of lower k. 

From Fig. [5] we note that at fixed values of the nor- 
malized soliton speed, Af , the amplitudes of the pertur- 
bations of the normalized potential, cool electron density 
and cool electron speed due to the solitary waves all in- 
crease with decreasing k. This is related to the increase of 
the true Mach number M/Mi for smaller k, as the phase 
velocity Mi is decreased. Thus larger disturbances are 
likely to be associated with increased suprathermality. 

Finally, turning to Fig. [HI two effects are observed: At 
fixed true Mach number, M/Mi, the soliton amplitude 
decreases with decreasing n (increasing suprathermality) . 
Despite that, the maximum values of soliton amplitude 
is found to occur not for a Maxwellian, but for the rel- 
atively low-K values of around 2.5-3. Thus, again, large 
observed amplitudes are likely to be associated with a 
low-K plasma. 

Hence, as shown above, these results could assist in 



the understanding of solitary waves observed in two- 
temperature space plasmas, which are often character- 
ized by a suprathermal electron distribution. 
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